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1. INTRODUCTION 

The neutron lifetime is still one of the least accurately measured fundamental nuclear constants. It is 
known that precision measurement of the neutron lifetime along with other neutron decay experiments 
has significant impact on the Standard Model parameters. Recent reviews are published in [1-3]. 

There are two principal ways to measure the neutron lifetime. In a neutron beam experiment the 
^-activity of a beam of cold neutrons is measured, and the neutron lifetime precision is determined by 
the accuracy of absolute measurement of neutron beam density in a well-defined fiducial volume in the 
beam and by the accuracy of absolute measurement of the neutron decay products (protons or electrons) 
count rate. 

In the ultracold neutron (UCN) storage method the neutron lifetime is measured by counting surviv- 
ing neutrons in a magnetic or material trap as a function of time. No absolute measurement is necessary 
in this case but serious problem arises in view of UCN losses in traps. In the UCN magnetic storage 
experiments the problem is in possible neutron spin- flip followed by its escape from the magnetic stor- 
age volume. Another problem is the marginal neutron trapping due to specular reflection of superbarrier 
neutrons from the ideally smooth magnetic mirror leading to appearance of quasi-bound neutron orbits. 

In material traps the main problem is to account properly for the UCN losses in their collisions 
with the surface of a material confinement cavity. One tries to minimize these losses choosing the 
materials with the lowest neutron capture and suppressing the neutron upscattering by lowering the trap 
temperature. 

The value of necessary corrections and systematic errors in inferring the neutron lifetime from the 
UCN storage data depend on these losses. 

The present experimental data base for the free neutron decay is large but contains significant discrep- 
ancy between data. The Table shows the values of the neutron lifetime obtained in the beam experiments 
and in the UCN storage experiments. Only the results with uncertainties not exceeding 10 s were taken 
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into consideration. 



Results of the neutron lifetime measurements in the beam experiments and in the UCN storage 

experiments 



Beam experiments 



Storage experiments 



891±9 (1988, P. E. Spivak) [4] 
*893.6±3.8±3.7 (1990, J. Byrne et al.) [5] 
889.2±3.0±3.8 (1996, J. Byrne et al.) [6] 
*886.8±1.2±3.2 (2003, M. S. Dewey et al.) [7] 
*886.6±1.2±3.2 (2004, J. S. Nico et al.) [8] 
886.3±1.2±3.2 (2005, M. S. Dewey et al.) [9] 



877±10 (1989, F. Anton et al.) [10] 
870±8 (1989, A. G. Kharitonov et al.) [11] 
887.6±3.0 (1989, W. Mampe et al.) [12] 
888.4±3.3 (1992, V. Nesvizhevsky et al.) [13] 
882.6±2.7 (1993, W. Mampe et al.) [14] 
885.4±0.9±0.4 (2000, S. Arzumanov et al.) [15] 
881.±3.0 (2000, A. Pichlmaier et al.) [16] 
878.5±0.7±0.3 (2004, A. Serebrov et al.) [17] 
* 874.61^6 (2007, V. Ezhov et al.) [18] 
878.2±1.9 (2009, V.Ezhov et al.) [19] 



Averaged value 


Average value without [17] and [19] 


887.6±2.7 


885.±0.82 




Average value including [17] and [19] 




881.3±0.53 



The data marked with the * are not used for obtaining the average values, only the last results of the 
corresponding experiment are taken into consideration. 

The world average value [20] based on [4,6,9,12,13,14,15] is 885.7±0.8 s. 

It is seen that the most precise UCN storage results may be divided in two groups: one - Refs.[12- 
15] and another - Refs.[17,19], with the work [16] between these groups. In view of this disagreement 
between UCN storage measurements and some difference between the results of the beam and the UCN 
storage experiments, new precision measurements of the neutron lifetime by different methods are de- 
sirable. It is essential that in the UCN storage experiments the total neutron disappearance probabiUty is 
measured, whereas the beam experiments are only sensitive to the neutron beta decay. 

2. DIFFERENT FORMS OF THE UCN MATERIAL TRAPS 

Different geometries (Fig. 1) of storage volume are considered here: cyUndrical with vertically and 
horizontally directed axis of the cylinder in respect to the Earth gravitational field, spherical, paraboloidal. 
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Figure 1 : The forms of the UCN traps. 
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conical, and vase. The quasi-spherical and horizontal cylinder geometries have been used in the neutron 
lifetime measurements [11,13,17], the vertical cylinder has been used in the first neutron lifetime experi- 
ments with UCN [21], the conical, vase and paraboloidal forms have never been used and are considered 
for comparison. 

In the presence of gravity, for equilibrium distribution of ultracold neutrons in the available phase 
space in a cavity, the neutron loss rate (s~^) is determined by the equation (we follow here the approach 
of [22]): 

d j Q{r, V, t)dh d\^- j //(v)(vns)^(r, v, t) d\ dS dt. (1) 

Here, Q{r, v, t) is the UCN phase space density, ng is the normal to the surface at the point of 
neutron collision with the wall. The neutron loss probability per one collision for isotropic neutron 
incidence, according to the accepted formalism is: 

p(v) = 2 77'^, K(y) = arcsin y - yJl-y'^, y = v/vb = Je/V, (2) 
y ^ 

where v and E are the neutron velocity and energy, respectively, Vb is the neutron boundary velocity of 
the cavity, rj is the UCN loss coefficient determined by the properties of the wall material: 

77 = - Im [// Re C/, U ^ {h'^/2m)47r ^ Nik, V = ReU, Imb = -a/2X. (3) 

i 

Here the complex potential U describes UCN interaction with walls, m is the neutron mass, Ni is the 
number density of nuclei of type i in the material, bi is the corresponding coherent scattering length on 
a bound nucleus of the wall, and a is the cross section of inelastic processes (capture plus upscattering) 
for neutrons with wavelength A. 

The neutron density depends on the height coordinate z defined relative to the lowest point of a 
cavity: 

g{z,v)^c5{vl-2gz-v^), (4) 

where c is the normalization constant, vq and ^(0, v) = c S{vl — v'^) are the neutron velocity and the 
neutron density at the bottom of the trap, g = 9.80665 m/s^ is gravitational acceleration. 
Solution of Eq. (1) gives the UCN loss rate for different forms of the traps. 

General expression for the UCN loss probabiUty (s~^) of the neutron with velocity v at the bottom of 
the cavity is 

where S{y) is the surface loss term and V{y) is the effective volume. 

In a cylindrical trap with horizontal axis the expressions for S{y) and V{y) are: 

S(y) = 5eyi(y) + ^SM, ^Rh£' KUy''-2gR{l-cos^)/vi\ d^, 
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^pi(y) = /J" K(^^Jy''-2gz/vi^ V2Rz - dz, 
V{y) = 2H £ sjy^ - 2gz/vl ^2Rz - z'^ dz, (6) 



where Scy\{y) is the loss contribution of the cylinder surface, Sp\{y) is the loss contribution of each of 
two plane surfaces, = arccos(^l — ^^^), Zq = {yvb)^/2g, R and H are the radius and the length of 
the horizontal cylinder, respectively. 

For a spherical trap of radius R these expressions have the form: 

S(y) = ttR^ K^^Jy^-2gz/vl^ dz, V{y) = £ ^y^ - 2gz/vl {2Rz - z^) dz. (7) 

For a conical cavity with the angle Q and flat bottom of radius ro the loss rate is determined by: 

S{y) = 5con(y) + ^bot(Z/), ^con(Z/) = VT ^ K y^ - 2gz/vi^ 

■77" r*^ fZQ I 

S^oi{y) = V{y) ^nj^ .Jy'' - 2gzlvl{z tan^)^ dz, (8) 

where ScoiXv) is the conical surface losses and 5'bot(y) is the bottom surface losses. 
For a vertical cylinder trap with the radius R the loss rate is determined by: 

S{v) = Scy\{y) + S^ot{y), Scyiiy) =t^ R i^(^y'^ - 2gz/vl^ dz, 

Sy.M^^K{y), V{y)=nR'J^ ^y' - 2gz/vl dz, (9) 



where Scy\{y) is the contribution of the cylinder surface. 

For a paraboloidal cavity (z — r^/R) the loss rate is determined by 

S(y) f\(^^Jy-'-2gz/vi^ ^Rz + dz, V(y) = nR _ 2gz/vl z dz. (10) 

For a cavity in the form of a vase {z — {RrY^'" — b) with the flat bottom of radius tq and radius 
at the height H the surface loss and the effective volume are: 



^va.e(y) = J £ 1^ Uy^ - 2gz/vi\ {z + 6)" 



S{y) = 5'vasc(y) + 5'bot(|/), 

n2(;2 + p)2("-i) 



sUv) = -Y<y)^ 



IT /"-^O I 

viy) = j^2l i^+pf'^^/y'-^gz/vidz. (u) 

Here, R = H''/{rU'' - rl^"") andp = (i?ro)^/". 
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Figure 2: The reduced UCN loss probability (s^ ) as a function of neutron energy for different forms of 
neutron traps: 1 - vertical cylinder, R = 50 cm, 2 - horizontal cylinder with radius of 50 cm and length 
100 cm, i - sphere, R = 50 cm, 4 - paraboloid z = r^/25, 5 -truncated cone, = tt/ 6, the bottom radius 
ro = 5 cm, 6 - truncated cone, 9 = 7r/4, the bottom radius ro = 5 cm, 7 - vase cavity {z = \fRr — p), 
with bottom radius of 5 cm, top radius of 50 cm and height of 100 cm, 8 - infinite plane. 
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Fig. 2 shows the reduced (the loss coefficient rj = 1) UCN loss probability (s^ ) as a function of 
neutron energy at the bottom of the traps of different form: cylindrical with vertically and horizontally di- 
rected axis of the cylinder in respect to the Earth gravitational field, spherical, paraboloidal, conical, and 
vase. For comparison the loss probability is shown also for neutrons bouncing at the infinite horizontal 
plane. 

AH results are shown for Vf, = 447 cm/s - the boundary velocity for the fully fluorinated poly ether oils 
(perfluoropolyethers - PFPE) used in the recent UCN lifetime measurement and having good neutron 
reflection properties [12,16,17]. 

As is seen, there is significant difference in the values and energy dependence of the neutron losses 
in the traps of similar size but different form. The traps expanding to the top, especially conical ones and 
in the form of vase, have low energy dependence of losses in a wide energy range. 

3. EXTRACTION OF THE NEUTRON DECAY CONSTANT 

The total loss rate r^^ of neutrons stored in a trap is determined by the neutron decay constant r^cay 
and the wall loss rate r^^^g. 

'^st — '^decay + '^loss" (1-^) 

There are two methods of extracting the neutron decay constant from storage measurements in the 
presence of losses. 

The first one consists in exclusion of the UCN surface losses from the experimental r^^ using several 
traps with the same UCN loss coefficient rj but with different UCN free path lengths between collisions 
with walls: 

= "^.Tecay + Vli, (13) 

where r^j and 7^ are the inverse experimental storage time and the calculated reduced (without rj) loss 
in i-th trap, respectively. 

The first method is realized introducing additional surface at the same storage volume [21] or using 
the traps of same form but different size (size extrapolation) [12,14,16,17]. 

For extrapolation the measured r^] are plotted against calculated 7^ or against calculated inverse free 
path length between neutron collisions with walls to obtain the intercept with the 7 = 0- axis. 

In the second method the UCN storage is performed in the same trap but with a sequence of different 
neutron energy spectra - the energy extrapolation method [12,13,17]. 

Sharp energy dependence of losses is essential for extracting the neutron decay constant from the 
storage data by the energy extrapolation method. Weak energy dependence of losses on the neutron 
energy gives (for example in the energy range from 30 to 100 neV for the curve 7 at Fig. 2) relative 
insensitivity of the storage time to the neutron spectrum. Conical and vase traps with their weak energy 
dependence of losses can be used more efficiently with the size extrapolation. 

One key systematic effect, that may arise when using the size extrapolation method with two or more 
neutron traps, is possible difference in the value of neutron losses at the surface of different traps. On 
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Figure 3: The scheme of the UCN neutron Ufetime measurement. 

the other hand in the energy extrapolation method, applied to the UCN storage data obtained in the same 
neutron trap, the homogeneity of the neutron reflecting coverage of the trap surface is not less important. 

Contrary to the size extrapolation the energy extrapolation needs more precise information on the 
realistic energy dependence of collision neutron losses. The energy dependence of UCN losses has been 
measured twice: first for copper surface in [23] and then with better precision for the Fomblin oil and 
grease surfaces in [24] . No deviation has been found in both measurements from the energy dependence 
of Eq. (2) except for the increase in losses for the neutrons with energy in the vicinity of the boundary 
potential of Fomblin oil [24]. This could be caused by additional small UCN upscattering by liquid 
surface modes of Fomblin. 

Actually, it is better to combine both methods of extrapolation in the neutron lifetime experiment in 
material traps [17]. 

The scheme considered here in more details is shown in Fig. 3. The storage measurements are 
performed in two storage traps of different size (e.g., two vertical cylinders with the diameters of 50 
and 100 cm) for a sequence of energy spectra of stored neutrons: with the upper energy cut-off from 
the lowest possible, say 10 neV, up to the largest possible for the used material, e.g. ~100 neV for 



8 



fluoropolymers. At small (100-200 s) storage time the initial integral spectra of stored neutrons is 
measured, at long storage (~ 1000-2000 s) the measured storage time comprises the neutron decay as 
well as the energy dependent losses. 

4. THE UCN SPECTRA FORMATION 

Sharp cutoff in the spectrum may be obtained with the help of a neutron absorber placed horizontally 
in the upper part of the storage volume in the cavity. The quality of the UCN absorber is determined 
by the height of the potential of the absorbing material. Traditionally, polyethylene (CH2)„ with its 
comparatively low negative potential for neutrons C/(CH2)n ~ — 9neV, was considered as an "ideal" 
absorber. 

Metal absorbers with low or negative scattering length like Ti or V may be not the best choice as the 
most effective UCN absorbers in view of oxidation of its surface. Even very thin - dozens A- oxide film 
forming positive potential barrier entails significant UCN reflection from the surface. Fig. 4 shows the 
neutron energy dependence of reflectivity at isotropic neutron incidence for polyethylene together with 
two other possible materials with good compensation of positive and negative coherent scattering lengths 
of nuclei - glyceride of melyssinic acid CggHisgOe with the height of the potential U = — 3.3neV and 
Zirconium hydride ZrH2 with almost perfect compensation of the scattering lengths: U = —0.325 neV. 

The characteristic "cleaning time" - the loss time for UCN with the energies exceeding the cutoff 
energy determined by the height of the neutron absorber, may be calculated in the same formalism as 
was used for calculations the neutron loss rate in traps (Eqs. (l)-(5)). 

The absorption probabiHty of the neutron with the energy E incident on a surface with an ideal step 
function profile at the angle 9 in respect to the surface normal is 

4 cos e ./cos2 e - (EJE) 
T{E, e, E,) = \ (14) 

(C0S^ + ^C0S2^- 

where £'5 is the boundary energy of the absorber. 

The expression for the absorption time constant for a vertical cylinder trap with the absorber at the 
height H with respect to the bottom of the cylinder, has the view: 

Kg 

Tclean(^^,-H") = V^, (15) 

where 

r./2Asm{29)cose{cos''9Tvl/{y'-2gH)f'\ 



COS 9 + (cos2 9 T vll{v^ - 2gH)) 



J<3 = 5 , (16) 

vl — 2\Ei,\/m, and the signs "-" and "+" correspond to positive and negative boundary energies of the 
absorber respectively. In the case of positive boundary energy of the absorber the integration over angle 
is performed over angles when cos^(^) > v^/ {v^ — 2gH). The Eqs. (16) are valid when > 2gH. 
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Figure 4: The energy dependence of the UCN reflection probability from polyethylene (CH2)n, glyceride 
of melissynic acid CgeHiggOe and Zirconium hydride ZrH2 at isotropic neutron incidence. 
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Figure 5: The energy dependence of the UCN cleaning time in a vertical cylinder with an absorber at 
H=50 cm for polyethylene, ZrH2 and an ideal absorber. The latter practically coincides with an ideal 
absorber. 

Fig. 5 shows the energy dependence of the "cleaning time" for three UCN absorbers: polyethylene, 
ZrH2, and an ideal absorber (100 % loss probability for all neutron energies) in the vertical cylinder with 
the absorber placed at the height H = 50 cm. The "cleaning times" for ZrH2 and an ideal absorber are 
practically the same. 

It is seen that the polyethylene absorber is generally not seriously inferior in respect to an ideal one 
except for the energy range close to the cutoff energy, where the cleaning time for the polyethylene may 
exceed several times the cleaning time for an ideal absorber. 

The time needed for cleaning spectrum should not be too large in view of additional neutron losses 
during cleaning process. 

In [25] it was proposed to construct a new gravitational trap similar to described in [17] - rotating 
horizontal cylinder with the radius of 120 cm. It is obvious that the spectrum cleaning time for such a 
trap may be very large. 

The UCN escape from a horizontal cylinder having the longitudinal hole with the angular width 26^0 
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Figure 6: Energy dependence of the calculated UCN escape time from horizontal cylinders: 1 - the 
diameter of 50 cm with the hole in the upper side with 29o = 0.4 rad, 2 - the half-cylinder of radius 50 
cm, 3 - the half-cyUnder of radius 120 cm. 

in the upper side of the cylinder may be imitated by placing the flat vertical ideal absorber as an upward 
directed prolongation of the hole. The escape rate is calculated according: 

f = T7T~^n7T^^ W = 7, / [v -2gz)dz, 

Vl{v) + V2{v) 2 JR(l+coseo) 

rn — do I 

Vi{v) =2HR^ / iAn^d^v'^ -2gR{l-coi^d)de, 

V^iy) = 2HR / _ 2gz dz, (17) 

Ji?(i+coseo) 

where S{v) is the loss term due to collisions with the absorber, Vi{v) and V2{v) are effective volumes of 
the assumed to be non-absorbing in this case cylinder and the absorbing prolongation, respectively. 

The neutron escape time from horizontal cylinders is shown in Fig. 6 for three cases: 29=0.4 rad, the 
cylinder radius 25 cm, and 26^0 = vr (open half-cylinder), the cylinder radius 50 cm and 120 cm. 

5. THE SIMULATION OF THE EXPERIMENTS 
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Figure 7: The results of the simulation of the neutron lifetime experiment in the vertical cylinder trap 
with radius 50 cm. The simulation of the neutron storage was performed for the Maxwellian spectra with 
the energy cut-offs at 20, 40, 60, 80, 90 and 100 neV. The neutron decay lifetime was 885.7 s. 

The results of simulation of the experiment with the UCN energy extrapolation is illustrated in Fig. 
7 for the vertical cylinder trap of the radius 50 cm. 

The storage times and the values of collision losses are calculated for 6 different cut-offs of the 
initial stored UCN spectrum: from 20 neV to 100 neV. The loss coefficient was taken r] = 2 x 10^® 
[26,17], the neutron decay lifetime 885.7 s [20], the neutron loss function has the form of Eq. (2), and 
it was assumed that initial neutron spectrum at each measurement is the Maxwell-like tail with perfect 
cut-off at the corresponding maximum energy. It was assumed that the storage times were measured 
with a statistical uncertainty of 4 x 10"'' (100 and 1000 s measurements, each with 10^ counts), and 
were distributed randomly. To reach this statistical precision needs ~ 530 cycles of fillings the neutron 
storage cavity if the number of stored neutrons in full spectrum with the i?rnax=100 neV is 10® per filling 
and linearly depends on the cutoff energy. 

In reality extrapolation of storage data to zero losses takes place in condition of non-complete infor- 
mation about the energy spectrum. The incident integral UCN spectrum is measured at small storage 
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Figure 8: The reduced UCN loss probability in the vertical cylinder trap for different UCN spectra as 
a function of reduced UCN loss probability for the Maxwellian tail spectrum in an assumption of the 
standard reflection law of Eq. (2). 

time by changing the height of the absorber. But precision of this measurement is not perfect and the 
incident UCN spectrum in the trap may differ from the Maxwell-like form. Therefore extrapolation to 
zero losses in an assumption of of the Maxwell-like spectrum may incur errors. 

Fig. 8 shows the calculated loss rate for different UCN spectra as a function of the calculated loss 
rate for the Maxwell-like spectrum. 

The extrapolated losses to zero energy differ from zero in the cases, when the spectra are not 
Maxwellian, by the value ~ 10^^ s^^ ^ 10^'^r^7^ if the r] = 2 x 10^^. It follows from this prelimi- 
nary estimate that the measurement of the incident UCN spectra with 5% precision will decrease this 
correction to the level 5 x 10^^ r^^. 

More serious correction may appear if the real dependence of UCN losses on the neutron energy 
differs from the form of Eq. (2) calculated for the ideal step potential. Analysis shows that the UCN loss 
energy dependence should be determined in this experiment with the precision no less than 5%. 

Fig. 9 shows the calculated UCN loss rate for the cases when there is 5% admixture to the standard 
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Figure 9: The same but for the Maxwellian tail spectrum and different UCN dependent UCN losses. 

energy dependence of Eq. (2), of different energy dependent terms in the UCN loss, as a function of the 
calculated UCN loss for the standard, Eq. (2) formula. 

The UCN spectrum was assumed to correspond to the tail of the Maxwell-like spectrum in these 
calculations. The extrapolated losses to zero energy differ from zero in the cases, when the neutron 
losses are not standard, by the value ~ lO^^s^^ ^ lO^'^ r^^ if the r] = 2 x 10^^ in the worst case of the 
square root UCN energy dependence of the UCN loss coefficient. 

It is obvious that for higher precision the loss corrections have to be evaluated more precisely and 
should be made as low as possible. 

The low loss probability per one neutron collision with the wall is expected from the result of recent 
experiment [17], where it was obtained r] ^ 2x 10^^ for the low temperature fluorinated oil (PFPE), and 
there is hope to decrease further this value below 10~^ [27,28]. On the other hand, larger size storage 
traps compared to [17] should decrease the loss rate due to a decrease of the wall collision rate. 

6. THE "ACCORDION" TRAP 

The effects of gravity on the storage of UCN in material traps with neutron losses have been consid- 
ered first in [22] and in [29] with implications for experiments to determine the neutron lifetime. The 



15 



consideration in [29] was focused on rectangular traps similar to the used in the experiments [12,16]. In 
these experiments the extrapolation to zero losses was performed by changing the length of the trap in 
horizontal direction. 

An interesting proposal in the paper [29] was to use the UCN trap in the form of a bellows with 
horizontal axis. When the length of bellows and consequently the volume of the trap is changed, the 
surface is remained unchanged. In this case extrapolation to zero losses is straightforward. The project 
to realize this approach is published in [30]. 

The model calculations of the UCN losses in wall collisions as a function of the length of the bellows 
are presented below. 

In the absence of gravity the accordion geometry is ideal in the sense of extrapolations to the infinite 
volume. In the presence of gravity the upper and lower parts of the trap are not equivalent, and the 
change of the form of the trap may, in principle, give some nonzero systematic effect for the value of 
extrapolated UCN losses. 

The exact geometric form of the surface for different stretches of the bellows is not known. Here two 
models are used: the linear model (truncated cones) and the sinusoidal one. One segment for the linear 
model of the bellows is shown in Fig. 10. 

It is assumed that at the stretching of the bellows the medium radius of the bellows is unchanged but 
the outer -Rmax and the inner i?min radii are changed according to the following expressions: 

d^H/N, sine ^d/l, r = i?o + -^(^- J), 

Z = r cos (/?, Rmax,min = -Rq ± (1^) 

2 tan 9 

where H is the length of the bellows, is the number of truncated conic segments of the bellows, d is 
the length of the segment in horizontal direction, Rq is the medium radius of the bellows, / is the length 
of the generatrix of the conic segment, 9 is the deformation angle, (p is the angle around horizontal axis 
of the bellows. In this case the datum in z direction is at the axis of the bellows. 
Similar to the calculation for other UCN trap forms we have: 

S{y) = N-S,on{y) + 2-SM, 

Scon{y) = [ dx [ ruljy'^ - 2gr cos ^p/vl] dip, 
sm 9 Jo Jo \ J 

SAy) = C^'^i^'Jy^-'^dz/vi^ ^j2Roz - dz, 

V{y) = 2N dx j"^ ^Jy^ - 2gz/vl \/2rz - z^ dz, (19) 

where Scon{y) is the surface loss on one conic section of the bellows, Spi{y) is the surface losses at the 
two disc surfaces of the bellows, V{y) is the effective volume term. 
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Figure 10: Schematic view of the "accordion" segment. 
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Figure 1 1 : Reduced UCN loss probability in the "accordion" trap for different UCN energies as a func- 
tion of the inverse volume of the trap. The datum is placed at the axis of the trap. 

The results of calculations of the loss rate as a function of the inverse volume of the bellows is 
shown in Fig. 11. The parameters of the trap were taken from the Ref. [30]: / = 2.9 cm, Rq = 
{Rraax + R-min)/'^ = 26.75 cm, the length was varied between 30 and 120 cm. 

Linear extrapolation 7(1/1^) = 7inf + &x {1/V) for each of the energies gave 7inf ~ (— 5±5) x 10"^. 
At the loss coefficient rj = 10^^ the systematic effect arising from the deformation of the UCN storage 
volume at the size extrapolation is about 5 x 10^^ s^^ for the inverse lifetime or ^ 0.04s for the neutron 
lifetime and does not depend on energy. 

At larger value of / the systematic effect at the extrapolation to infinite volume is increased. 

The serious problem in the use of bellows UCN trap is deposition of the low temperature coverage at 
the corrugated surface and its stability at the deformation of the bellows. Another disadvantage consists 
in the narrow energy range of confined neutrons ^ 30 — 70 neV for the datum placed at the axis of the 
trap and in the case of the PFPE coverage for the bellows size of [30]. These limits follow from the 
conditions that the UCN should have the energy sufficient to reach the highest points of the trap and not 
exceed the boundary energy of the PFPE Eb ^ 106 neV. 
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7. MAGNETIC TRAPS 

One of the most interesting methods of slow neutron storage has been proposed and experimentally 
demonstrated by Paul and collaborators [31,10]. It was a toroidal hexapole magnetic trap for very slow 
neutrons in the /xeV energy range: the neutrons with longitudinal velocities between 8 and 15 m/s and 
radial velocity spread up to 4 m/s were stored in the ring with the diameter of ~ 110 cm. The neutron 
injection from the neutron guide into the magnetic ring was performed with the help of fast pneumatically 
driven neutron mirror reflector. The final experiments, in which the neutron lifetime was measured with 
the precision 10 s, were performed at the very cold neutron channel of the PF-2 neutron source [32] of 
the ILL High Flux Reactor. 

The great advantage may be achieved using pulsed injection of neutrons from the pulse neutron 
source [33,34]. In this method the neutron bunch from the UCN converter after spreading over the 
neutron guide can be trapped in the UCN cavity with high efficiency. 

The detailed MCNP calculations for a solid para-hydrogen converter with an optimal cold ortho- 
hydrogen premoderator located in the radial channel of the MARK-2 TRIGA reactor yield the production 
of ~ 10^ neutrons/Mole/ 10 MJ reactor pulse in the neutron velocity half-cylinder with the radius of 4 
m/s and the longitudinal velocity interval from 8 to 15 m/s . Possible size of one-mole para-hydrogen 
converter: the diameter of 7 cm, the thickness of 0.5 cm is optimal for the reactor radial channel and for 
the neutron extraction in this energy range. 

At the distance L from the neutron source to the entrance to the storage ring, and the neutron velocity 
interval of stored neutrons between vi and f 2, the length of the neutron bunch passed this distance is 
L{{v2/vi) — 1). The length of the neutron guide at the TRIGA reactor from the moderator to the entrance 
of the magnetic ring may be as short as about 5-6 m/s, and at the length of the storage ring ~350 cm 
[31] practically optimal condition may be realized for the neutron injection: almost all the bunch in the 
longitudinal velocity range 8-15 m/s fills the storage ring during one reactor pulse. With larger storage 
ring and consequently larger orbit time the injection and storage will be more effective. 

The neutron injection efficiency is not reported in [31]. In an assumption that the neutron survival 
during transport through the neutron guide and injection is only 1% it is expected up to 10^ neutrons 
trapped in the ring per one filling. It is three orders of the magnitude larger than reported in [31]. 

The additional advantage consists in the possibility to variate the spectrum of injected neutrons with 
help of the pulse injector due to time-of-flight separation of neutrons in the neutron guide. 

Similar possibility of pulsed neutron filling from pulse neutron source exists for other types of the 
magnetic traps for ultracold neutrons [35] if one provides the magnetic trap with pulse magnetic shutter. 

It seems possible, in principle, to control neutron losses from magnetic ring by placing very slow 
neutron detector based on the low temperature scintillator (for example pure Csl [36]) at the outer side 
of the ring storage volume and PMT in the central part of the ring. 

Another possibility - the decay proton counting is possible placing proton counters in the center of 
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the storage ring and introducing accelerating electric field between the walls of the storage volume and 
the proton detectors. 

The author is grateful to Dr. V. K. Ignatovich for clarifications concerning his work [22]. 
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